Nonlinear, ground-state, pump-probe spectroscopy 
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A theory of pump-probe spectroscopy is developed in which optical fields drive two-quantum, 
Raman-like transitions between ground state sublevels. Three fields are incident on an ensemble of 
atoms. Two of the fields act as the pump field for the two-quantum transitions. The absorption or 
gain of an additional probe field is monitored as a function of its detuning from one of the fields which 
constitutes the pump field. Although the probe absorption spectrum displays features common to 
those found in pump-probe spectroscopy of single-quantum transitions, new interference effects are 
found to modify the spectrum. Many of these features can be explained within the context of a 
dressed atom picture. 
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I. INTRODUCTION 



Of fundamental interest in nonlinear spectroscopy is 
the response of an atomic vapor to the simultaneous ap- 
plication of a pump and a probe field. A calculation of the 
probe field absorption is relatively straightforward [QJ|] 
in the weak probe field limit. Let fl and £1' denote the 
pump and probe field frequencies, A = Q — lu the pump 
field detuning from atomic resonance w, and Si = Q' — O 
the probe-pump detuning. For a pump field detuning 
|A| 3> 7 e , Xj where 7 e is the upper state decay rate and 
X is a pump-field Rabi frequency, one finds the spectrum 
to consist of three components. There is an absorption 
peak centered near <5i = —A (O' = to), an emission peak 
centered near Si — A (O' = 2S1 — ui) and a dispersive like 
structure centered near Si = 0. Experimentally, a spec- 
trum exhibiting all these features was first obtained by 
Wu et al. H . The absorption and emission peaks can be 
given a simple interpretation in a dressed-atom picture 
j|, but the non-secular structure centered at Si = is 
somewhat more difficult to interpret ^,^|. The width of 
these spectral components is on the order of j e , neglect- 
ing any Doppler broadening. 

The spectral response can change dramatically when 
atomic recoil accompanying the absorption or emission of 
radiation becomes a factor M, as in the case of a highly 
collimated atomic beam or for atoms cooled below the 
recoil limit. In this limit, the absorption and emission 
peaks are each replaced by an absorption-emission dou- 
blet, and the dispersive- like structure is replaced by a pair 
of absorption-emission doublets. The spectrum can be 
given a simple interpretation in terms of a dressed atom 
theory, including quantization of the atoms' center-of- 
mass motion |7j. It turns out, however, that at most one 
absorption-emission doublet (one of the central ones) can 
be resolved unless the excited state decay rate is smaller 
than the recoil shift. Since this condition is violated for 
allowed electronic transitions, it is of some interest to 
look for alternative level schemes in which this structure 



can be resolved fully. If the optical transitions are re- 
placed by two-photon, Raman-like transitions between 
ground state levels, the widths of the various spectral 
components are determined by ground state relaxation 
rates, rather than excited state decay rates. As a result, 
the probe's spectral response should be fully resolvable. 
Raman processes have taken on added importance in sub- 
Doppler H and sub- recoil [ | cooling, atom focusing [jl0| , 
atom interferomet ry pd| - [l4 |, and as a method for probing 
Bose condensates |I5|], 

In this article we propose a scheme for pump-probe 
spectroscopy of an atomic vapor using Raman transi- 
tions. This is but one of a class of interactions that can be 
considered under the general heading of nonlinear ground 
state spectroscopy. The spectral response is found to be 
similar to that of traditional pump-probe spectroscopy 
|EJ; however, new interference phenomena can modify the 
spectrum [Sec. III]. The interference phenomena can be 
interpreted in terms of a dressed atom picture [Sec. IV]. 
Although part of the motivation for this work is the study 
of recoil effects, such effects arc neglected in this article. 



II. EQUATIONS OF MOTION 



The atom field geometry is indicated schematically in 
Fig. |l|. Three-level atoms interact with two optical fields, 
Ei and E2, producing strong coupling between initial 
and final levels I and 2 via an intermediate excited state 
level e. Field Ei couples only levels 1 and e, while field 
E2 couples only levels 2 and e. In addition, there is 
a weak probe field E that couples only levels I and e. 
As a consequence, fields E and E2 can also drive two- 
photon transitions between levels 1 and 2. Levels 1 and 
2 are pumped incoherently at rates Ai and A2, respec- 
tively, and both states decay at rate T. The incoherent 
pumping and decay represent an oversimplified model for 
atoms entering and leaving the interaction volume. 
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FIG. 1. Schematic di- 
agram of the atom-field 
system. Fields Ei and E 
drive only the 1 — e tran- 
sition and field E2 only 
the 2 — e transition. 



The incident fields are assumed to be nearly copropa- 
gating so that all two-photon Doppler shifts can be ne- 
glected. In this limit and in the limit of large detuning 
on each single photon transition, one can consider the 
atoms to be stationary with regards to their interaction 
with the external fields. We wish to calculate the lin- 
ear probe absorption spectrum. The electric field can be 
written as 



E(R,t) = - 



£ ie *(ki-R-nit) + £ 2e i(k 2 -R-r22t) + £. e i(k-R-nt) 



+ c.c. 



(1) 



where fii, f2 2 , and f2 are the field frequencies, ki, k2, 
and k the field propagation vectors, and c.c. stands for 
complex conjugate. In an interaction representation, ne- 
glecting any decay or incoherent pumping of the ground 
state levels, the state probability amplitudes obey the 
equations of motion. 



-i (7e/2) a e , 



idi = xie tAlt a e 



ia 2 = X2e 



(2a) 
(2b) 
(2c) 



where Xj = —d e jEj/2h (j = 1, 2) and x = —d e iE/2h are 
Rabi frequencies (assumed to be real and positive), d e j is 
a dipole moment matrix element, and Aj = Qj — u> e j and 
A = il — ui e i are atom-field detunings. Assuming that 
the magnitude of the detunings are much larger than 7 e 
and any Doppler shifts associated with the single photon 
transitions, it is possible to adiabatically eliminate the 
excited state amplitude to arrive at the following equa- 
tions for the ground state amplitudes: 



»6i = S1O1 + S (e lSlt + e~ lSlt ) 01 + ge~ l5t a 2 
id 2 = S 2 a 2 + ge lSt ai + g'e l5>t ai, 



where 



5 = A 2 - Ax = fi 2 -Oi +W21; 
6' = A 2 -A = n 2 -n + w 2 i; 
Si = A- A x = 0-Qi = 5-6', 



(3a) 
(3b) 

(4a) 
(4b) 
(4c) 



are detunings associated with two-quantum processes 
and 

g = X1X2/A; g' = XX1I A; 
S 1 =Xi/A; S 2 =x\IA- S = X Xi/\ (5) 

are Rabi frequencies or Stark shifts associated with two 
quantum processes. In writing Eqs. (ph, we assumed that 



A w Ai w A 2 and |A| > 



It will prove convenient, especially when going over to 
a dressed atom picture, to introduce a representation in 
which 



Oi - 6ie -«*/2 e -i(S 1+ S 2 )t/2. 

Combining Eqs. (||) and (^) one finds 



(6) 



ib x = -{5/2)b x + gb 2 + S (e lSlt + e~ lSlt ) b x 

+g'e l5lt b 2 ; (7a) 
^i 2 = (5/2)b 2 +gb 1 +g'e- Ult b 1 , (7b) 



where 



5 = 6-(S 1 -S 2 ). 



(8) 



The corresponding equations for density matrix elements 
Pu = \bi\ 2 , P22 = \b 2 \ 2 , pu = b x b* 2 = p* 21 are 



Pu = ~ig (P2i - P12) ~ ig'e l5lt p 2 i + ig'e- tSlt p 12 - Tp u + A l5 (9a) 
P22 = ig {P21 - P12) + ig'e i5lt p 21 - ig'e~ Ult p 12 - Tp 22 + A 2 ; (9b) 
P12 = i6pi 2 - ig ( P22 - pu) - ig'e Ult (p 22 - Pll ) - iS (e l5lt + e"^ 1 *) p 12 - Tp 12 , (9c) 



where the incoherent pumping and decay terms have 
been introduced. It is important to note that, in this 
representation, the frequency appearing in the g' terms 
is 61 = 5' — 6 = O — Oi . In other words, the effective field 
frequency associated with field E 2 in this representation 



is f^i rather than Q 2 . 

It follows from the Maxwell-Bloch equations that the 
probe absorption coefficient, a, and index change, An, 
are given by 
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where N is the atomic density, 



Pie 



1 

A 



(0) 

XPn 



XiPu 



(10a) 
(10b) 

(11) 



and pf^ , , and p\ 2 are coefficients that appear in the 
solution of Eqs. (||) (to first order in \) written in the 
form: 



(0) 



3,3 



1,2 (12) 



The first and third terms in Eq. (]Tl]) are analogous to 
the terms that appear in conventional theories of pump- 
probe spectroscopy, but the second term is new and leads 
to qualitatively new features in the probe absorption 
spectrum. 

An expression for p' le is given in Appendix A. The ab- 
sorption coefficient is plotted in Figs. ||(a)-(c) for several 
values of 8/g, and 



n 



= \Jxi/Xi 



(13) 



If T) -C 1, the two-quantum probe absorption spectrum 
has the same structure as the probe absorption spec- 
trum involving single quantum transitions. The situa- 
tion changes if r\ > 1. For example, aside from an inter- 
change of absorption and gain components as a function 
of 8, the probe spectrum for single quantum transitions 
depends only on the magnitude of the pump field detun- 
ing. This is clearly not the case for two-quantum tran- 
sitions, as is evident from Fig. ||(a) drawn for rj = 1, 
T/g = 0.1, 8/g — ±1. Probe absorption and gain are 
interchanged when 8 changes sign, but the ratio of the 
amplitude of the absorption to gain peak changes when 8 
changes sign. There is another subtle difference present 
in these spectra. The sense of the central dispersive com- 
ponent is opposite to that for single quantum transitions. 
With decreasing 77 , the sense of the central component 
would reverse, as the spectrum reverts to the same struc- 
ture found in pump-probe spectroscopy of single quan- 
tum transitions. The probe response also depends on 
the sign of A (through g = X1X2/A); this feature follows 
from the dependence of the spectrum on the sign of 8 and 
the relationship 



Pi 



,(-8,-A,-S 1 ) = ~p' le (S,A,S 1 y 



(14) 



which can be derived using Eqs. (A3)-(A7) of Appendix 
A. It is also possible for the components centered at pos- 
itive or negative to vanish (in the secular approxima- 
tion) for certain values of 77, as can be seen in Fig. 0(b). 
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(b) 



-1.5 



0.0 
s/g 



8/a=0;I7q=0.1m=1/5l 




l8/a=0; r/q=0.1 ; n=1 1 



-3.0 -1.5 0.0 1.5 3.0 



-3.0 -1.5 0.0 1.5 3.0 
s,/g 




(C) 
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FIG. 2. Probe field absorption in arbitrary units. Positive 
ordinate values correspond to probe absorption and negative 
values to probe gain. 

The case of 8/g = is shown in Fig. |^(c) for r\ = 
1/5,1,5, and A > 0. If 77 = 1/5, the spectrum is simi- 
lar to that found for single quantum transitions M. For 
77 = 1, the spectral component at negative Si is found to 
vanish. When 77 > 1, there is a dispersive-like structure 
centered at 81 = that is not found in the pump-probe 
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FIG. 3. Dressed-state energy level diagram. In the interac- 
tion representation adopted in the text, the frequency of field 
E2 must be set equal to fii in calculating resonance condi- 
tions. For (Aa — As) > 0, solid arrows correspond to probe 
absorption centered at Si = luba and dashed arrows corre- 
spond to probe gain centered at 81 = —loba- 

It is now possible to use the energy level diagram (Fig. 
||) to read directly the probe absorption spectrum. The 
probe field is absorbed (or amplified) via two quantum 
transitions between states \A) and \B). The two quan- 
tum transitions involve one photon from the probe field 
and one photon from either field E\ or E2, since all of 
these fields couple states \A) and \B) to state \e). It is 
important to remember that the effective field frequency 
of field E2 is equal to f^i in this interaction representa- 
tion. Fields Ei and E couple state |e) to the components 
of states I A) and \B) involving state |1), while field E 2 
couples state |e) to the components of states \A) and \B) 
involving state |2). For example the matrix element for 
the two-quantum process from state \A) to \B) involving 
absorption of a probe photon and emission of a field E2 
photon is 



-IX 
-iA 



cos(0) 



(-%) 



r - i{5\ - lo B a 



■ cos(0), 



while that for absorption of a probe photon and emission 
of a field Ex photon is 



cos(0) 



Hxi) 



r - i(Si - luba 



■ipsin(O). 



These two processes add coherently, such that probe ab- 
sorption via transitions from state \A) to \B) is propor- 
tional to the sum of these two matrix elements squared, 
multiplied by the population difference [p^a ~ Pbb) • ^ n 
other words, the probe absorption at <5i = luba is pro- 
portional to a quantity C+ given by 



C+ = (g/AT) [(Aj - A 2 ) /T] cos (26) ( ^ sin (0) cos(0) + - cos 2 (0) 

V 



(21) 



Similarly, probe gain via transitions from state \A) to \B) at 5\ = —luba is proportional to 



4 



C_ = (g/AT) [(Ai -A a )/T]cos(20)) ( iprj sin (0) cos (0) - - sin' (0) 



(22) 



A formal derivation of these results is given in Ap- 
pendix B. 

For the sake of definiteness, let us take (Ai — A 2 ) > 0; 
then C+ corresponds to absorption for 5 > and to gain 
for 5 < 0, while C_ corresponds to gain for 5 > and to 
absorption for S < 0. Note that the component centered 
at Si — —lob a vanishes if A > and tan(0) = r/ 2 , while 
that at Si = luba vanishes if A < and tan(0) = r\~ 2 . 
The values of A± = ±C*± [r 2 A/ \g\ (Ai - A 2 )] are plot- 
ted in Fig. | as a function of 5/g for A > and 
r) = 1,2. For A < 0, one can use the relationship 
A±(-A,-«)=A F (A,«). 
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FIG. 4. Amplitude A+ of the peak centered at Si = uiba 
and amplitude A- of the peak centered at Si = — cjba, for 
A > 0. Positive values of A± correspond to absorption and 
negative values to gain. 



The probe absorption vanishes in the secular approx- 
imation ( |l9| ) when 5 = 0, since, in this case, 6 = j 
and the populations of the dressed states are equal. The 
lowest order dressed atom approach is not useful in this 
limit. Typical spectra are shown in Fig. [2j(c) and were 
discussed in Sec. III. 



IV. CONCLUSION 



The probe absorption spectrum has been calculated for 
two-quantum transitions between levels that are simulta- 
neously driven by a two-quantum pump field of arbitrary 
intensity. In addition to features found in conventional 
pump-probe spectroscopy of single quantum transitions, 
new features have been found that can be identified with 
interference phenomena. Both Doppler and recoil effects 
were neglected in out treatment. For nearly copropagat- 
ing fields, effects arising from these processes are negligi- 
ble. Doppler shifts can be accounted for by the replace- 
ments Si — » Si + (ki — k) ■ v, Si — S — » Si— 5+ (k 2 — k) • v, 
and S\ + 5 — » 8\ + 5 + (2ki — k 2 — k) • v in the equations 
in the Appendix. 

The dependence of the interference effect of the signs 
of A and S can be understood in the bare atom picture 
in a perturbative limit. A schematic representation of 
the probability amplitude leading to probe absorption at 
5\ = 8 is shown in Fig. ^(a). Each arrow represents an 
interaction with one of the fields. The two contributions 
to the final state amplitude add coherently. Putting in 
the appropriate energy denominators, one finds that the 
absorption varies as 



A = 



i 2 XX*2 



■4 * I |2 

i XX2 \Xi\ 



( 7e /2-iA)[r-i(A-A 2 )] ( 7e /2-iA)[r-i(A-A 1 )]( 7e /2-iA)[r-i(A-A 2 ) 



(23) 



For 



^> r, and |A| ^> j e , this equation reduces to 

9 



.4 



XX 2 



A 



1 



r 2 + (<Ji 



~\2 

5 



1 



\9\Tfi> 



(24) 



which shows the dependence on the signs of A (tp = 
| A | / A) and S. A similar calculation for the emission com- 



ponent represented schematically in Fig. ||(b) leads to 



G 



* 2 * 

X X1X2 



A 2 5 



r 2 



(25) 



New effects will arise if the fields are not copropagating 
and the active medium is a subrecoil cooled atomic 
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FIG. 5. Schematic representation of the 1— > 2 tran- 
sition probability leading to probe absorption or probe 
gain in lowest order perturbation theory in the bare ba- 
sis. The thin arrow represents the probe field, the broad 
filled arrows field E\, and the broad open arrows field Ei- 
(a) absorption, (b) gain. Terms involving the sequential 
absorption and emission of the same field have been ne- 
glected, since such terms result only in Stark shifts of lev- 
els 1 and 2. The diagrams are drawn for 8 > 0; if S < 0, 
the roles of absorption and gain would be interchanged. 



vapor, a highly collimated atomic beam, or a BEC. As 
for single quantum transitions , each component of the 
spectrum undergoes recoil splitting. Since the center-of- 
mass momentum states differ for two-quantum processes 
involving fields E\ and E from those involving fields E% 
and E, one might expect the spectrum consists of eight 
absorption and eight emission components rather then 
the four absorption and four emission components found 
for single quantum transitions; however this does not ap- 
pear to be the case. Instead, each component results 
from a coherent superposition of two quantum processes 



involving fields {E\,E) and (E2,E). 
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APPENDIX A: BARE STATE CALCULATIONS 

Substituting Eqs. ( [l2| ) into Eqs. @, one finds to ze- 
roth order in the probe field that 



T 2 + S 2 



r 2 + S 2 + A\g\ 



(o) 
Pxi 



(A 2 + Ai) 



,(o) 



(o) -ig (o) / jo) 

912 r - is ~ V 21 



and that, to first order in the probe field, w~* 
Pi2' <°2i' an d m+ = P22 + P11 sa tisfy 



0: 



(r + iS{) w+ - 2ig*p+ 1 + 2igp+ = 2ig'p$- 



[r + i(5i- 5)] p+ 2 + ig*w + = -ig' 
[r + i(5i + 5)] p^\ — igw + = iSp$ . 

Equation (O) can be rewritten as 



iSpl 



P'u « (X7A) 



(0) 1 
P11 ~ -jVo + V12 



jj (Ala) 

(Alb) 

(Ale) 

P22 ~ P11 1 

(A2a) 
(A2b) 
W> ; (A2c) 
(A2d) 

(A3) 



where 



Va = (xi/xT w + ; 
Vi2 = (Xa/xT 

,' X2X* \ + 
V21 = — P21- 

X XiJ 



(A4a) 
(A4b) 

(A4c) 



In Eqs. @-@) we have allowed the Rabi frequencies to 
be complex. 

The quantities ?/o, 2/12, and 2/21 satisfy the coupled 
equations: 

[r + i{5x- &)] 2/12 + i \g\ tpr]- 2 y = a; (A5a) 

(r + %Si)yo - 2i \g\ Wv2i + 2i \g\ Wvxi = b\ (A5b) 

[r + * (<5i + 5)} y 21 - i \g\ V^J/o = c > ( A 5 C ) 



where 



b = zlkUm = 2c 

r + is 

a=-i \g\^~ 2 w^ -JiL-wW. 

r — is 



(A6a) 
(A6b) 



and ip = (|A| / A). Note that the equations do not de- 
pend on the phase of the various Rabi frequencies, but 
do depend on the sign of A. Explicit solutions for 2/0 and 
J/12 are: 



G 



Vo = -« 



2a \g\ ipT] 2 (5 + St - iT) + b 


5 2 - - iT) 2 


+ 2c \g\ i>rf (6-81. + iT) 


(5 1 - iT) (s 2 -5 2 + 2iT5 l + T2 + 4 \g\ 2 ^ 





5 2 +5 [Sx - iT) - 2ir<5i - T 2 -2 \g\ 2 - br/- 2 \g\ tp (6 + Si - iT) - 2c \g 



Vi2 = i- 



The line shape is totally non-secular when 6 = 0. In 
the limit that A > 0, \g\ > T, and r <C r] 2 , one finds 
that the absorption coefficient a for 8i ~ is 

1 ( kNd 2 e \ Aj-A 2 5jT 

^-iI^aJ r <*; + !*) ' (A8a) 

that for 6i w 2 |g| is 



(5i - ir) (,5 2 - <5? + 2zr5 x + r 2 + 4 |. g | 2 ) 

Xi = IXi|e # \X2 = lx 2 |e^,x=lx|e^. 



(A7a) 
(A7b) 

(B5) 



1 ( kNd 2 le \ Ai-Ag (fr-2| g |)r 



\2he Q AJ T 



and that for <$i ss —2 |g| is 



(5 1 -2\g\f+T 2 



1 f kNd 2 e \ Ax-A 2 (ji + 2jgj)T 



8 \2he A J T 



(6 1 + 2\g\) 2 +T 2 



(1 + T 2 ), 



(A8b) 



(l- V - 2 ). 



(A8c) 



Note that the component at 5i = —2 \g\ vanishes if rj = 1. 
For A < 0, one can use Eq. (111). 



APPENDIX B: DRESSED-STATE 
CALCULATIONS 

Equation (0) can be written in the form 
ihb = (V + V p ) b, 

where 

-5/2 g* 



v = n 



Vn = h 



9 5/2)^ 
Se i5lt + S*e~ lSlt g'*e lSlt 



g' e - Ult 







_ X1X2 1 _ XX2 Q _ X Xi 

9- 9 -£-> 

and we have allowed for complex Rabi frequencies, 



(Bl) 



(B2) 



(B3) 



(B4) 



If one introduces semi-classical dressed states via the 
transformation 



where 



T r b, 



cos (0) e 



i<t>d/2 



ai<Pd/2 



-i<Pd/2 s 
i<Pd/2 . 



sin (0) 



and 



sin (0) er^l 2 cos (0) 



62+2(1-^) 



(B6) 



(B7) 



(B8) 



(B9) 



(recall that -0 = |A| / A), then the dressed-state Hamil- 
tonian is given by 



v d = h 



-cjsa/2 



ujba/2 
The dressed state density matrix, 



T c V p Tt. 



PAA PAS 

Pba Pbb 



evolves as 



A A 
A B 



(BIO) 



(Bll) 



(B12) 



Off-diagonal terms have been neglected in the matrix rep- 
resenting the incoherent pumping, since they give rise to 
terms of order T/u>ba <C 1 (secular approximation). 
The dressed state density matrix is expanded as 



P d = pf+P>^ 



(B13) 



and it is found from Eqs. that pj 

obeys the equation of motion 



d_ 

dt 





(ujba + 5i)pba 



{uJBA-5i)p\ B 





1 

ift 



(B14) 



where 



7 



(9) [Scoa {9) - g'*e 1 ^ sin (9)\ ; cos (0) 
in (0) [-g'* sin (0) + S cos (0)] ; sin (0) 



cos 
sin 



Ssin(0) + g'* cos (0) e^ d ] 
Ssin (0) + g'* cos (6)e^A 



In the secular approximation, the steady state solution of Eq. ( B14| ) is 

pi- 



PIb 
Pea 



where 



p+ B = i cos (0) [S sin (9) + g>* cos (9) e*'] (pf A - ) / (T + i - WflA )) , 
p+ A = -t sin (9) [-g'* sin (0) e 4 *' + 5 cos (0)] (p^ - p^) / (T + i (fc + WflA )) • 



The coherence p' le needed in Eq. (10) for the absorption coefficient and index change is given by 



Pie « (1/A) 



* (0) . * + , * + 

X Pii + XiPn + X2P12 



(B15) 

(B16) 

(B17a) 
(B17b) 

(B18) 



The first term can be evaluated using Eq. (Alb) for p^; it contributes to the index change, but not the absorption. 
For the remaining terms, one rewrites and p^ 2 in the dressed basis using Eqs. (B6),flB"g|),(Bll), and uses Eq. (|B5|) 
to extract all the phase factors to arrive at 



p' le «(x7|A|) Mi } +/+ + /■ 



where 



/+ = 
/- = 



*\9 



[r + i (St - lu ba )\ 
i\g\ 



cos (29) 



(Al A2) cos 2 (9) ( ipr] sin (9) + - cos (6) 



r 



■ cos (26») 



^ Al sin 2 {9) ( ijm cos (9) - - sin (9) 



[v + i {St + lu ba )\ "~ y -"' r 

Note that the approach and results of Sec. Ill are unchanged if one uses complex dressed states defined by 



1^4) 
\B) 



|1) 

|2) 



(B19) 

(B20a) 
(B20b) 

(B21) 
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